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Abstract 


The  bandstructures  of  the  layer  compounds  lT-TiS2  and  lT-TiSe^  have 
been  calculated  within  the  local  density  approximation  using  the 
self-consistent  LAPW  method.  Both  compounds  are  found  to  be  semimetals 
with  T  -L  overlaps  of  0.24  and  0.55  eV,  respectively.  The  positioning 
and  occupation  of  the  bands  near  the  Fermi  level  are  not  consistent  with 
the  carrier  densities  determined  from  Hall  measurements;  however 
agreement  is  obtained  when  the  overlap  is  decreased  by  about  0.25  eV. 
Calculations  performed  with  contracted  lattice  constants  reproduce 
qualitatively  the  same  trends  in  Hall  coefficients  as  observed 
experimentally.  However,  quantitative  agreement  can  not  be  obtained 
without  an  overlap  decrease  of  about  0.25  eV  -  indicating  that  TiSe? 
is  a  semimetal  and  TiS2  a  semiconductor. 

*Prnsent  address:  Physics  Department,  3aylor  University,  Waco,  Texas  76793 

U.  S.  A. 
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I .  Introduction 

The  IT  layer  compounds  Ti$2  and  TiSe^  have  been  intensively  investiga¬ 
ted  for  nore  than  a  decade,  not  only  because  of  their  technological  uses  in 
high  energy-density  batteries,  but  also  because  of  their  nearly  or  partially 
overlapping  conduction  and  valence  bands.  The  question  of  their  intrinsic 
seniconducting  or  semimetall ic  bandst ructures  has  presented  a  challenge  to 
both  experimental  and  theoretical  techniques.  The  original  optical 
exper  irr.ent  s  on  TiS^  and  TiSe?  suggested  that  both  were  serai  conductors  with 
fundamental  gaps  of  1-2  eV  (Greenaway  and  Nitsche  1965);  the  metallic 
properties  of  each  being  attributed  to  inpurities  and  defects.  Tsbeuchi 
and  Katsuta  (1970a,  1970b)  proposed  that  the  metallic  properties  were 
intrinsic,  resulting  from  the  overlap  of  the  chalcogen  p  and  titanium  d 

bands.  For  TiS2  this  view  has  been  supported  by  extensive  resistivity 
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experiments  on  highly  pure  samples  which  have  revealed  a  T  dependence  over 
a  range  of  10-400K.  (Thompson  1975);  by  studies  of  the  variation  of  the 
conductivity,  therraopower ,  and  lattice  parameter  with  stoichiometry 
(Thompson  et.  al^  1972);  and  by  angle-integrated  photoenission  experiments 
which  have  shown  a  long  tail  in  the  density  of  states  indicating  band 
overlap  or  a  band  gap  of  no  more  than  0.1  eV  (Williams  and  Shepherd  1973, 
Fischer  1973).  Angl e-resolved  photoenission  (AKF)  experiments  have  been 
more  ambiguous  in  that,  while  they  have  shown  overlapping  bands  in  TiSe2  and 
a  small  gap  in  TiS2,  only  a  few  bands  at  high-syrametry  points  have  been 
resolved  (Bachrach  £t  al_  1976,  Traum  et_  a_l_  1978,  Chen  e_t_  al^  1980,  Barry  e£ 
al  1983). 


More  recently,  Mall  coefficient  measurements  and  reflectivity  data, 

taken  together,  have  established  a  strong  case  for  a  seniconcucting  picture 

for  TiS^  and  a  seraimetallic  one  for  TiSe^  (Friend  -_et_  al_  1977,  Logothetis  e_t 

al  1979,  Kukkonen  £t_  al_  1981);  while  further  resistivity  measurements, 

2 

showing  a  deviation  from  the  previously  reported  T“  behaviour,  have  eroded 

much  of  the  TiS^  semimetallic  support  (Kukkonen  et_  a_l_  1981,  Klipstein  et  al 

1983)-  However,  within  the  TiS^  serai  conductor  model,  the  origin  of  about 
20  3 

10  electons  per  cm  in  the  highest  purity  samples,  remains  unknown. 
Wilson  (1978)  has  suggested  that  excess  Ti  atoms  might  contribute  more  than 
four  electrons  to  fill  the  bands,  but  so  far  this  has  not  been  supported  by 
any  evidence. 

To  try  to  clarify  the  experimental  and  theoretical  pictures  for  these 
compounds,  we  have  calculated  sel f-consistently  the  bandstructures  of  TiS0 
and  TiSe^  using  the  Linearized  Augmented  Plane  Wave  (LAPW)  method.  A  self- 
consistent  (SC)  approach  is  essential  in  examining  these  compounds  since  the 
early  non-SC  calculations  gave  gaps  that  were  far  too  large:  2.0  eV  TiS2, 
3.5  eV  TiSe^  (Murray  and  Yoffe  1972);  1.4  eV  TiS2>  0.5  eV  TiSe2  (Myron  and 
Freeman  1974).  More  recent  SC  calculations  have  come  closer  to  experiment: 
0.23  eV  gap-TiS2,  0.18  eV  overlap-TiSe2  (Zunger  and  Freeman  1977,  1978);  0.5 
eV  overlap-TiS2  (Umrigar  et^  1982);  0.S0  eV  overlap-TiS2 ,  0.95  eV 

overlap-TiSe2  (Temroerman  et_  a_l  1983).  In  the  present  work,  we  find 
senimetallic  bandstructures  for  both  TiS2  and  TiSe2>  with  F  to  L  band 
overlaps  of  0.24  and  0.55  eV,  respectively. 

We  have  also  repeated  the  bandstructure  calculations  using  contracted 
lattice  constants  in  order  to  understand  the  differing  behaviour  of  the  TiS2 
and  TiSe2  Hall  coefficients  under  pressure.  We  find  the  TiS2  Hall 
coefficient  to  be  much  less  pressure-dependent  than  TiSe2,  thus 
qualitatively  reproducing  the  experimental  trend.  However,  in  order  to  get 


quantitative  agreement,  the  evidence  points  to  a  reduction  in  overlap  by 
about  0.25  eV  from  that  calculated. 


II.  Theoretical  Method 


The  LAPV.'  method  used  here  has  been  well  described  in  the  literature 
(Andersen  1975,  Koelling  and  Arbman  1975).  Briefly,  the  LAPW  basis 
functions  are  formed  by  dividing  real  space  into  an  interstitial  region  and 
muffin-tin  (MT)  spheres,  then  approximating  the  potential  and  solving  a 
one-electron  equation  in  each.  The  resulting  solutions  are  joined 
continuously  and  dif  f  c-rent  iably  at  the  region  boundaries.  Although  a  MT 
potential  is  used  to  construct  the  basis  functions,  so  long  as  the  basis  is 
sufficiently  complete  it  may  be  used  in  a  variational  calculation  employing 
a  potential  having  no  shape  approximations.  Inside  the  MT  spheres  a 
seni-relativistic  form  of  Dirac's  equation  is  solved  using  a  spherically 
averaged  potential  (Koelling  and  Harmon  1977).  Both  the  radial  solutions 
and  their  energy  derivatives  are  included  in  the  basis  to  enable  the  matrix 
elements  to  be  independent  of  the  energy  ever  the  range  of  interest.  In  the 

i 

interstitial  region  the  basis  functions  are  the  usual  plane  waves. 

The  total  Hamiltonian  H  can  be  written  as  the  sum, 

H  -  V  +  MT!  +  AVNS’ 

where  contains  the  kinetic,  exchange-correlation,  and  MT  potential 

terms,  AV^  is  the  correction  potential  due  to  non-MT  terms  in  the  intersti¬ 
tial  region,  and  AV^,<,  is  the  correction  potential  due  to  non-spher ical  terms 


inside  the  MT  spheres. 


In  the  present  work  we  have  included  AV^  exactly  and  A\\. ^  approximately 
using  a  scheme  proposed  by  Andersen  (Jepsen  et_  a_l  1982).  In  this 
approximation  the  plane  wave  expansion  of  the  potential  is  continued  inside 
the  muffin-tin  spheres  and  the  spherical  average  subtracted  from  it.  The 
new  plane  wave  expansion  differs  significantly  from  the  old  only  in  the 
region  just  inside  the  MT  sphere  boundary.  (Symmetry  causes  the  new 
expansion  to  be  small  near  the  lattice  sites.)  Since  the  LAPW  basis 
functions  are  continuous  and  have  a  continuous  first  derivative  across  the 
boundary,  the  interstitial  plane  waves  may  be  used  to  evaluate  the 
(approximate)  non-spherical  matrix  elements.  In  testing  this  approximation 
for  0  ,  where  the  non-spherical  corrections  to  a  varped-MT  potential  are 
sizable,  Umrigar  found  that  about  90«  of  the  eigenvalue  error  was 
eliminated.  Such  an  approximation  should  be  even  better  in  a  solid  where 
the  higher  coordination  number  reduces  the  deviation  from  spherical 
symmetry.  We  have  therefore  utilized  this  approximate  full-potential 
throughout  the  present  work. 

In  iterating  to  self-consistency  we  have  used  a  basis  set  consisting  of 
approximately  150  plane  waves.  Exchange  and  correlation  were  treated  using 
the  Kedin-Lundqvist  (1971)  potential.  For  our  Brillouin  zone  integrations 
in  each  iteration  we  have  used  a  set  of  twelve  special  k-points  for  D 

id 

symmetry  derived  from  the  six  Ic-point  set  for  symmetry.  For  the 
densities  of  status,  we  have  calculated  eigenvalues  over  a  uniform  mesh  of 
75  k-points  in  the  irreducible  wedge  of  the  Brillouin  zone.  Self-consistency 
was  considered  to  be  achieved  when  input  and  output  potentials  differed  by 
less  than  12  mRy;  the  eigenvalues  being  converged  to  better  than  2  mRy. 


111.  Results 


In  figures  1  and  2  are  displayed  our  calculated  bandst ruciures  of 
stoichiometric  TiS7  and  TiSe^.  along  with  the  correspond ing  densities  of 
states.  Both  compounds  are  found  to  be  seraimetals- with  T  to  L  band  overlaps 
of  0.2A  and  0.55  eV ,  respectively.  In  addition  there  is  an  indirect  F-M 
overlap  of  0.3A  eV  in  TiSe7,  but  no  similar  overlap  in  TiS^-  The  Fermi 
energy  in  HS2  falls  just  above  a  deep  minimum  in  the  DOS,  in  agreement  with 
photoemission  spectra  (Williams  and  Shepherd  1973). 

The  strongly  overlapping  (0.55  eV)  conduction  and  valence  bands  in 

TiSe^  represent  the  first  theoretical  agreement  with  the  overlap  observed  in 
the  photoemission  experiments  of  Bachrach  et_  al_  (1976).  However,  although 
the  Ti  d-band  was  not  always  resolved  unambiguously,  two  other  ARP  studies 

of  TiSe.,  report  a  smaller  band  overlap  of  about  0.3  eV  (Traura  £^t_  jal^  1973, 

Chen  et  al  1980).  These  conflicting  experimental  results  have  not 
previously  been  reconciled.  If  only  the  latter  two  studies  are  accurate, 
then  our  calculated  conduction  band  in  TiSe^  is  too  low  by  about  0.25  eV. 
(We  shall  return  to  this  possibility  in  the  discussion  below.)  However,  a 
reconciliation  of  all  three  ARP  experiments  is  possible  within  our 

bandstructure:  both  the  F~  and  valence-band  maximum  r*  states  lie  near  the 
Fermi  level  and  are  separated  by  0.2  eV.  The  0.3  eV  overlap  may  then 
correspond  to  the  overlap  between  the  and  L^  states,  while  the 

measurement  of  Bachrach,  e_t_  al_  on  less-stoichiometric  samples  reveals  the 
full  F*  to  l|  semimetallic  band  overlap. 

The  dimensions  of  our  calculated  electron  pocket  at  L  (.22  Lll  x  .2S  LA 

>:  1,M)  arc  in  excellent  agreement  with  the  experimental  values  of  .2  LH  x  .3 

LA  (Traura  e£  a_l  1978)  and  .2  LH  x  .A  LA  (Chen  et_  al_  I960)  on  higher  purity 

TiScn  samples.  However,  none  of  these  pocket  sizes  are  consistent  with  the 
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carrier  concentration  of  about  2  x  10  /cm  determined  from  Hall 
mcasui  cncnt  s  on  the  best  Tif^  samples  (Friend  et^  al_  19S2).  Our  calculated 


L  pocket,  corresponding  to  a  concentration  of  6.6  x  10  /era  ,  is  thus  at 

least  a  factor  of  three  too  large.  If  the  band  overlap  is  decreased  by  6.1:3 

eV,  an  L  pocket  of  size  .15LK  x  .19LA  x  .70LM  results,  corresponding  to  2.2 
20  3 

x  10  carriers  per  era  .  This  is  further  evidence  that  the  overlap  cay  be 
overest iraated  by  about  0.25  eV. 

A  similar  overestiraation  of  the  band  overlap  cay  occur  in  our 

calculation  for  TiS2.  As  previously  stated,  we  find  the  bandst ructure  to  be 

seninetallic  with  an  overlap  of  0.24  eV,  in  contradiction  to  the  strong 

exper icental  semiconductor  support.  The  positioning  of  the  valence  bands 

agrees  well  with  recent  angle-resolved  photoeraission  experiments  (see  Figure 

4  of  Bairy  et  al  1983).  Our  calculated  electron  pocket  size  of  .15  LH  x  .19 

LA  x  .47  LM  is  in  fair  agreement  with  Barry  _et^  _al' ' s  highest-purity  sample 

pocket  size  of  .06  LH  x  .17  LA,  since  the  ARP  results  depend  upon  the 

contour  criterion  chosen  (this  at  half -maximum  intensity)  and  the  angular 

acceptance  of  the  analyser  (2°  corresponding  to  .06  LH) .  Unlike  TiSe^,  the 
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carrier  concentration  corresponding  to  this  occupation  (1.6  x  10  )  is  very 

20 

close  to  that  determined  from  the  Hall  coefficient  (1.4  x  10  ).  However, 

if  TiS2  is  truly  a  semiconductor ,  then  the  origin  of  these  carriers  in  the 
experimental  samples  must  be  impurities  or  structural  defects  -  not  the  band 
overlap  found  here. 

Earlier  TiS2  ARP  measurements  found  an  electron  pocket  of  size  .2  LH  x 
.4  LA  (Chen  et_  al_  1980);  but  it  was  supposed  that  this  larger  electron 
occupation  was  also  due  to  the  presence  of  impurities  and  defects.  If  we 
increase  the  occupation  of  our  electron  pocket  at  L  by  raising  the  Fermi 
energy  (corresponding  to  the  addition  of  impurity  or  defect  atoms),  then  we 
find  that  the  Fermi  level  is  pushed  entirely  out  of  the  valence  band  before 
a  pocket  of  size  equal  to  the  ARP  data  is  reached.  Evidence,  then,  of  only 
one  carrier  type  in  the  less-stoichiometric  TiS_  samples  is  not_  definitive 


in  characterizing  the  bandstructure  as  semiconducting  (Lcgcthetis  et_  al 
1979,  KuV.kcnen  et  zl_  1931).  Hc'-ever,  similar  evidence  on  much  purer  samples 
does  provide  evidence  that  TiS2  is  either  a  semiconductor  or  a  small  band 
overlap  (less  than  0.15  eV)  seminetal. 

In  tables  1  and  2  we  compare  the  present  results  with  those  of  previous 
calculations.  Our  measured  conduction  and  valence  band  widths  are  in  good 
agreement  with  those  of  Zunger  and  Freeman  (1977,  1978)  using  a  self-consis¬ 
tent  LCAO  method.  Our  calculated  direct  and  indirect  gaps  and  overlaps  are 
intermediate  between  Zunger  and  Freeman’s  and  those  calculated  by  Temmerman, 
et  al(13S3)  using  the  self-consistent  L>irO  method.  Both  our  results  and  those 
of  Temmerman  et  ai  are  semi-relativistic,  while  Zunger  ana  Freeman's  are  r.cn- 
relativistic . 

It  is  interesting  to  note  that  although  the  majority  of  early  calculations 
have  tended  to  over-estimate  the  valence  band  (VB)  -  conduction  band(C' 
aration  (seme  by  as  much  as  1  eV) ,  the  most  recent  self-consistar.t  calculations 
have  tended  to  swing  the  other  way.  The  gap  now  appears  to  be  underestinx *.ed 
not  only  in  TiS„  and  TiSe?,  but  in  other  complex  compounds  such  as  ?c?  (’Zimmer 
et  al  1979)  and  the  rare-earth  mor.cpnicrides  (Hcsegawa  1980).  It  is  in  r.o  way 
clear  that  the  continued  disagreement  with  experiment  as  to  the  exact  measure 
of  band  gap  or  overlap  is  due  to  the  inferiority  of  present  methods.  Instead, 
the  problem  may  be  related  to  the  well-known  underestimation  of  semiconductor 
band  gaps  (by  40,'i  or  more)  within  the  local  density  approximation  (see  Ferdew 
and  Levy  1983,  and  references  therein). 

In  order  to  understand  the  effects  of  pressure  upon  the  cne-electron 
bandstructure,  we  have  performed  calculations  using  contracted  lattice 
constants  for  both  TiSn  and  TiSe?.  Since  the  ccmpressiblit ies  of  these  two 
ccr.nour.es  have  not  been  measured,  we  have  substituted  these  of  the  lavered 


compounds  NbS?  and  N’bSe^  (Jones  j2_t  al_  1972).  Mien  these  compress  ibilities  are 
extrapolated  to  20  kbar  ('.here  the  calculations  were  performed),  cr.e  attains 
contractions  of  .  3271  and  2.271  for  the  TiS?  a  and  c  lattice  spacings,  and  .1271 
and  3.271  for  TiSe?.  Although  we  cannot  expect  the  N'b  compounds  to  reproduce 
exactly  the  TiS^  and  TiSe„  compressibilities,  we  do  expect  the  S  to  Se  tr-_nd 
(less  compressible  S,  more  easily  ccmpresscd  Se)  to  be  closely  the  same. 


Our  contracted  bandstructures  are  displayed  with  dashed  lines  in 
figures  3  and  A.  To  illustrate  the  changes  due  to  increasing  pressure,  we 
have  shifted  the  energy  scales  so  that  the  lowest  conduction  bands  align  at 
F.  The  most  interesting  feature  of  these  bandstructures  is  the  fact  that 
there  is  almost  no  variation  either  in  the  p-h  overlap  or  in  the  shape  of 
the  lowest  conduction  band  in  TiS„,  but  that  both  a  modified  conduction  band 
and  an  increased  overlap  (of  0.2  eV)  are  observed  in  TiSe.,. 

The  nature  of  these  differences  may  be  understood  by  comparing  the  four 


bands  at 

F  nearest 

the  Fermi  energy. 

In  both 

vneenpret'  sec  cor.pouncs 

the 

lowest  ly 

ing  state 

is  the  which  :• 

'  v. :v'St/d 

■primarily  of  chalcogen 

P, 

orb ital s 

ccmb in  ed 

to  form  a  state  :l 

'*;t  is  Lr. 

tibending  both  within 

and 

between  layers.  Extending  into  the  i n t  •  ■  r  1  ..y r  region  as  it  does,  this  state 
is  the  nose  sensitive  to  any  reduction  in  the  c  lattice  spacing.  In  both 


TiS 

7  and 

TiSe0  the  main. 

effect  of  t h e 

2C  kb or  of  p 

ressure  is  the 

shift  of 

thi 

5  band 

relative  to  the  other  three. 

In  TiS.,,  the 

shift  is  just 

up  to  the 

3“ 

state 

at  the  top  of 

the  valence  b;: 

.nd .  However 

in  TiSe„  the 

2  shifts 

through  the  valence-band  maximum  3  state  to  become  the  top  valence  band  in 
the  ccmoressec  structure.  (The  ordering  of  these  hands  was  confirmed  by  setting 
the  LAPW  enerttv  parameters  equal  to  the  Fermi  energy.  Thus,  inside  the  Mi's 
the  semi-relativistic  equations  were  solved  at  precisely  the  energy  of  interest, 
as  in  the  A?V  method.)  It  is  primarily  this  shift  of  the  2  above  the  3+  state 


that  causes  the  increase  in  the  TiSe,,  band  overlap.  This  difference  ir.  2 
positioning  in  the  two  compounds  is  maintained  even  when  the  bands  overlap 
less  strongly  than  our  calculations  here  indicate.  (Vie  determined  this  by 
repeating  the  calculations  using  an  unphysical  Xu  exchange  potential  with 
c=  1 . ) 

Interpreted  in  terms  of  the  Hall  coefficient,  these  results  shew  that 
with  the  Fermi  energy  near  the  tep  of  the  valence  band  in  (impure)  TiS.,,  the 
pressure-induced  valence-band  shifts  have  little  effect  on  the  measured  Hall 
coefficient  (<5%) .  In  contrast,  because  the  TiSe„,  Femi  energy  lies  within 
the  valence  band  even  with  similar  numbers  of  impurities  or  defects, 
pressure-induced  shifts  increase  the  electron  occupation  at  L  causing  the 
Hall  coefficient  to  decrease  over  the  entire  pressure  range  of  20  kbar.  At 
higher  pressures  we  expect  the  TiS2  ^  state  to  push  through  the  Fermi  level 
so  that  the  carrier  concentration  will  begin  to  rise  and  the  magnitude  of 
the  Hall  coefficient  to  decrease  in  much  the  sane  way  as  observed  in  TiSe?. 
However,  for  nearly  pure  TiS.,  samples,  the  percentage  change  in  Hall 
coefficient  surpasses  101  -  in  disagreement  with  experiment  (Friend  et_  al^ 
1977).  Thus,  even  though  we  observe  a  much  smaller  pressure  variation  of 
the  band  occupancy  in  TiS^i  we  are  not  able  to  account  for  the  measured 
constancy  of  the  Hall  coefficient.  Thus,  -the  difference  in  the 
pressure-dependence  of  the  Hall  coefficient  between  TiS,,  and  TiSe^  cav  not 
entirely  be  attributed  to  differences  in  the  band  shifts  of  two  semimetals; 
the  evidence  points  to  T1S2  being  a  semiconductor  and  TiSe^  being  a 
semimetal . 

In  summary,  we  have  calculated  the  electronic  bandstruc tures  of  TiS^ 
and  TiSe?  self-consistently  within  the  local  density  approximation,  finding 
them  to  be  semimetals  in  agreement  with  other  recent  calculations  (t'mrigar 
et  al  1  982  and  Temmerman  e_t_  a_l  1983).  The  potential  has  been  treated 
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accurately,  and  ve  believe  that  any  residual  shape  approximation  has  little 
efrect  on  the  band  overlap.  The  positioning  and  occupation  of  the  bands 
near  the  Fermi  energy  are  not  consistent  with  the  carrier  densities 
determined  from  Kali  measurements,  however,  good  agreement  is  obtained  when 
the  band  overlap  is  decreased  by  0.25  eV.  Pressure-induced  shifts  in  the 
semimetallic  bandstructure  reproduce  qualitatively  the  sane  trends  as  the 
Hall  data;  however,  quantitative  agreement  cannot  be  achieved  without  a 
decrease  in  the  calculated  overlap.  This  analysis  points  to  TiS.,  being  a 
semiconductor  and  TiSe0  being  a  semimetal. 
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Table  1  : 

Comparison 

of  TiS2 

band  widths  and  gaps 

ll 

obtained  by  different  self-consistent  (SC) 

and  non-SC 

methods 
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Table  2:  Comparison  of  TiSe^  band  widths  and  gaps 
obtained  by  different  self-consi stent  (SC)  and  non-SC  methods 
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re  Cantions: 


1.  TiS9  bandstructure  and  total  density  of  states. 

2.  TiSe?  bandstructure  and  total  density  of  states. 

3.  TiS?  bandstructure  near  the  Fermi  energy  calculated  with 
uncompressed  (solid  line)  and  compressed  (dashed  line) 
lattice  constants.  Energy  divisions  are  in  electron  volts 

4.  TiSe„  bandstructure  near  the  Fermi  energy  calculated  with 
uncompressed  (solid  line)  and  compressed  (dashed  line) 
lattice  constants.  Energy  divisions  are  in  electron  volts 
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